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In this paper we found an Exact solution for massless scalar field with cosmolog- 
ical constant. This exact solution generalized the Levi-Civita vacuum solution^ to 
CN ■ a massless scalar field, with a cosmological constant term. This solution in the ab- 

CU sence of the Cosmological constant recovers the spacetime of a massless scalar field 

^ . with cylindrical symmetry (Buchdahl metric [2]). Also if the scalar field disappears, 

the spacetime is a representation of de-Sitter space. We prove that the form of the 
metric's function which was purposed in [l| is valid even if we assume a general 
(— i ■ form. Too we show that in which conditions this solution satisfies energy conditions. 

Q-e Finally the validity of focusing theorem is proved. 

+^ ■ PACS numbers: 04.30.-w, 96.10.+i, ll.25.-w 
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INTRODUCTION 



In the [l[ the authors acclaimed that they were succeed in solving the problem of find- 
ing an exact solution of Einstein field equations for a massless scalar field with cylindrical 
symmetry in the presence of a cosmological const ant [TJ. They obtained a new two parameter 
exact solution (LB meric) 1 which recovers at least one member of Buchdahl family ^ in 
the absence of Cosmological constant and also LCA family{3j].Theire work in [1] completed 
previous works as [3, 6, 8, 11]. In that work the authors solved the field equations in a special 
case and no thing was stated about the general exact solution. Now we applied a method 
for solving field equations and show that we can construct a general two parameters ex- 
act solution which in the absence of the cosmological constant recover Buchdahl solution. 
We show that this system of differential equations posses only those solutions that the au- 



thors stated in [lj. Later a class of solutions of Einstein field equations is investigated for 
a cylindrically symmetric spacetime when the source of gravitation is a perfect fluid [4] .As 
a historical note we added here that, our solution with Cosmological constant and scalar 
field is a generalization of Levi-Civita family [7, 8]. Previously the exact solution of Einstein 
field equations with a cosmological constant term were found by Linet [10] and in a more 
efficient form by Tian [11]. The singularity problem in a family of cylindrically symmetric 
spacetimes which is described by the collapsing of scalar fields was discussed by Wang and 
Frankel [9].Senovilla presented an explicit exact solution of Einstein's equations for an inho- 
mogeneous dust universe with cylindrical symmetry [5]. In this work first we prove that the 
general solution for field equations only is one, which was stated in [lj.Then we investigated 



This is an abbreviation to the memories of H. A. Buchdahl and the presence of the cosmological constant 
term 
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the validity of Energy conditions and focusing theorem for our solution. Also a short section 
is present which discusses the stability of exact solution under small perturbations via via 
Lyapunov exponents method[27,28,29]. 

LFIELD EQUATIONS 

We begin with a general cylindrically symmetric metric in Weyl coordinates (t,r,ip,z), 

ds 2 = -e u(r) dt 2 + dr 2 + e v ^dip 2 + e w(r) dz 2 (1) 

Field equation for a massless minimally coupled scalar field in the presence of a cosmological 
constant term A is reading as 2 : 

Rltv - A 9l^u = <P;t,<t>;u (2) 

We labeled metric functions as Ui = {u(r),v(r),w(r)},(j) = </>(r) , h — In terms of these 
functions we can rewrite field equation (2) in the following succinct forms: 

3 

2< + u' t Uj — 4A = 0,2 = {1, 2, 3} (3) 
j'=i 

3=1 3=1 

Now we write equation (3) in the following simple form: 

|>^) = 2Ae^ = {l,2,3} (5) 

E 3 - «■ 

Where in it / = — 3 ^ 3 . If we do summation on i = 1, 2, 3 in (5) we can write a differential 
equation for new function / = f(r), 

±(fef) = SAef (6) 
In terms of this new function the equation (4) converted to the following equation: 

3/" + f' 2 - 7 A = <p 12 (7) 
Substituting /" from (6) in (7) we obtain the following integral for 0: 

= ±V2 J dry/ A -ft (8) 

Thus if we solve the equation (6) we can obtain both metric functions and (f>. 



I will mostly use natural units h = c = 1 and 8irG = 1. 
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About the stability of the field equations 

Now we investigate the stability of the field equations (3,4) under small pertubations. First 
we note that instead of working with the system (3,4) we can treat the (6,8) as the field 
equations. in linear approximation (6) becomes 

/" « 3A 

The general solution for this simple ODE is 

f( r ) = ~Y r + ci 

If we take the de-Sitter radius as a = yj~j^ this function is nothing but the usual metric 

function of the Asymptotic de-Sitter. If we want to check the eigenvalues of the linearized 
matrix for system of field equations we must change the system (3,4) to a higher rank first 
order system. For this, we introduce the following set of new variables, 

Ui = xd = {1,2,3},0 = y 

Then we have 

3 

2x'i + Xi^2xj-AA = 0,i = {l, 2, 3} 

3=1 

2E4 + E*?-4A = 4,' 2 

3=1 3=1 

For a stationary point in the phase space we must set all first order derivatives equal to zero, 

3 

x i x i ~ 4A = °' i = t 1 ' 2 ' 3 ^ 

3=1 

3 

^^ 2 -4A = 0,y = 

3=1 

In the first equation if we summing on the indices and comparing both of them we obtain 

^ = 2/<M = {l,2,3},y = 

We perturb the field equations as 

xi = 2/a + 5xii = {1,2,3},?/ = 5y 

Expanding the field equation in these perturbations up to first order and taking to the mind 
to having an asymptotic stability the real part of the eigenvalues must be negative. In the 
case of our field equations the matrix posses negative purely real eigenvalues, thus according 
to the Lyapunov theorem this system is stable. 
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ILEXACT SOLUTIONS 

In this section we investigate all possible solutions for (6), (8). 
The exact solution for ordinary differential equation (6) is 3 : 

/(r) = -V3Ar + i log^^e 2 ^ - c 2 ) 2 ) (9) 

We keep the relation (9) as a constraint for all solutions of (6).Substituing this function in 
(5) we can write the following general form of metric functions: 

We note here that if the solutions (10) satisfy in (9) , then we have Y^j=i a j = O.Thus the 
final form of metric functions in : 



2 cie 2v ^ 3Ar — Co 
Ui = - log( ^ 2 ) +P»i= 1, 2, 3 (11) 

Later we will determine the coefficients C\,C2- The condition Y^j=i@j = ~ §l°s(12A) may 
be satisfied by adding a constant (3j to each metric function which can be absorb in the 
non radial coordinates. In the languages of potential theory we can choose a new gauge 
for potential functions «j(r) .The reason is that in cylindrical spacetimes we know that 
the solving of the Einstein field equations can be reduced to solving a system of potential 
equations [3]. A better option is the comparison of (11) with the previous metric function in 



cqi 

0= 



Ui = u(r) = ±^r + \ log(l + ee 2Var ),i = L 2, 3 (12) 

As it was attend that in [l[by the authors , only the minus sign in metric function is consistent 
with scalar field equation of motion □</> = .In terms of a new parameter £ = ~i\J~^ s £ 
which was defined in [lj]- we must choose the constant (3j as the following: 

& = ~log(-C2),i = 1,2,3 (13) 
Obviously if C\ G R, then c 2 G C (purely imaginary) and consequently /3j is real. 



IIIrVALIDITY OF ENERGY CONDITION IN LB METRIC 

The origin of the null energy condition ('NEC) and of the strong energy conditions (SEC) 
is the Raychaudhuri equation together with the requirement that the gravity is attractive 
for a spacetime manifold endowed with a metric 4 [24]. 
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logo*) = JT f 

This is a Tetrad's representation of the metric 
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The classical energy conditions of general relativity, to the extent that one believes that 
they are a useful guide [12, 13], allow one to deduce physical constraints on the behavior 
of matter fields in strong gravitational fields or cosmological geometries. These conditions 
can most easily be stated in terms of the components of the stress energy tensor T^ v in an 
orthonormal frame. Ultimately, however, constraints on the stress-energy are converted, 
via the Einstein equations, to constraints on the spacetime geometry. 5 For a perfect fluid 
cosmology, and in terms of pressure and density, the so-called Null, Weak, Strong and 
Dominant energy conditions reduce to [14]: 

NEC:pi + p > 

WEC:This specializes to the NEC plus p > 
SEC: This specializes to the NEC plus Yl Pi + P > 
DEC: p^\pi\. 

Note particularly that in FRW models of Universe, the condition SEC is indepen- 
dent of the space extrinsic curvature k. Now, DEC implies WEC implies NEC, and 
SEC implies NEC, but otherwise the NEC, WEC, SEC, and DEC are mathematically 
independent assumptions. In particular, the SEC does not imply the WEC. Violating the 
NEC implies violating the DEC, SEC, and WEC as well [14]. Note that ideal relativistic 
fluids satisfy the DEC, and certainly all the known forms of normal matter encountered 
in our solar system satisfy the DEC. With sufficiently strong self-intereactions relativistic 
fluids can be made to violate the SEC (and DEC); but classical relativistic fluids always 
seem to satisfy the NEC. Most classical fields (apart from non-minimally coupled scalars) 
satisfy the NEC. Violating the NEC seems to require either quantum physics (which is ' 
unlikely to be a major contributor to the overall cosmological evolution of the universe) or 
non-minimally coupled scalar fields (implying that one is effectively adopting some form of 
scalar-tensor gravity). Using this dynamical formulation of the energy conditions, Santos 
et al. [15] derive some bounds, for the special case k = 0(flat cosmology), on the luminosity 
distance oIl of supernovae, and then contrast this with the legacy [16, 17] and gold [18] 
datasets. In reference [19] bounds on the distance modulus are presented for general values 
of k while in reference [20] they concentrate on the lookback time. 

Due to the lack of satisfactory dark energy models, many model-independent methods were 
proposed to study the properties of dark energy and the geometry of the universe[21,22,23]. 
Another very interesting and model-independent approach is to consider the energy 
conditions [24,25]. Recently Energy Conditions and Stability in f(R) theories of gravity 
with non-minimal coupling to matter is discussed and determined the bounds from the 
energy conditions on a general f{R) functional form in the framework of metric variational 
approach. [26] . 

IV: CHECKING THE FOCUSING THEOREM 

In this section first we state the focusing theorem in the manner of congruences both 
in time like and null cases. Then we checked that our metric which is constructed from 



5 In particular in a FRW spacetime one is ultimately imposing conditions on the scale factor and its time 
derivatives 
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a messless cylindrically symmetric scalar field in the presence of a cosmological constant 
satisfies this theorem or not?. 

a:Focusing theorem in General Relativity 

Let a congruence of time like geodesies be hypersurface orthogonal, means there is exist 
a vector field u a (time like, space like or null and not necessarily geodesic). We called it a 
hypersurface orthogonal if Co> Q/ g = U[ a . t p] = O.It means that there exist a scalar field $ such 
that u a oc $ a and let the (SEC) hold. So that from Einstein field equations 

R aP u a uP ^ 0. 

From Raychaudhuri equation [8[ implies: 

^ = ~\e 2 - a^a a p - R aP u a u? ^ (14) 

Where in it, 9 = u a a = B% is the expansion scalar , a a p = B^ a ^ — \9h a p the shear tensor, 
h a/ 3 the transverse part of g a p [si (which is purely 'spatial'). The expansion must therefor 
decrease during the congruence's evaluation. Focusing theorem stated that an initially di- 
verging (9 > 0) congruence will diverge less rapidly in the future, while an initially converging 
(9 < 0) congruence will converge more rapidly in the future. The physical interpretation is 
that gravitation is an attractive force when the (SEC) holds and the geodesies get focused 
as a result of this attraction. 

For a congruence of null geodesies be hypersurface orthogonal and (NEC) hold, too we have 
^ ^ where 9 = k a a and k a is the tangent vector field. It is important to realize that 
the Raychaudhuri equation is purely geometric and independent of the gravity theory under 
consideration. The connection with the gravity theory comes from the fact that, in order to 
relate the expansion variation with the energy-momentum tensor, one needs the field equa- 
tions to obtain the Ricci tensor. Thus, through the combination of the field equations and 
the Raychaudhuri equation, one can set physical conditions for the energy-momentum ten- 
sor. The requirement that gravity is attractive imposes constraints on the energy-momentum 
tensors and establishes which ones are compatible. Of course, this requirement may not hold 
at all instances. Indeed, a repulsive interaction is what is needed to avoid singularities as 
well as to achieve inflationary conditions, and to account the observed accelerated expansion 
of the universe. 

b : Straightforward calculations 

Now, we consider a congruence of radial, marginally bound, time like geodesies of the LB 
metric 

ds 2 = dr 2 + w{r){-dt 2 + d V 2 + dz 2 ) (15) 

Where, 

w (r) = e ~ 2 ^ r (ee 2VMr + 1) 2/3 
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For radial geodesics,the components of_4- vector field u a , u v = u z = 0, and the geodesies are 
marginally bound if — = —u t = E. This means that the conserved energy is precisely 



equal to the rest-mass energy, and this gives us the equation u l = ^y. 6 From LB metric |l| 
we know that w(r) > 1. Indeed since the radial like component of 4-vector velocity must be 
a real function one can deduced that 

\E\ > 1 

From the normalization condition 

g a pu a uP = -1 

we have: 



//'■ ,1— 1 



w[r\ 



The upper sign applies to outgoing geodesies, and the lower sign applies to ingoing 
geodesies. The 4- velocity is given by: 




u a = (—-±\l^-- 1,0,0) (16) 
w[r) y w(r) 

And , using (16) we can write: 

U a d a = ^-a t ± (4^T - l) 1/2 ^r 

w[r) w[r) 
E 2 

u a dx a = -Edt ± (— - - l) 1/2 dr 
w(r) 

It follows that u a is equal to a gradient of a scalar function $ where: 

U a = a 

and: 

$ = _t± [fJ?-iy/2 dr (17) 

J w(r) 

The integral could be written in terms of hypergeometric functions which we don't write it 
here .This expression means that the congruence is everywhere orthogonal to the spacelike 
hypersurfaces $ = constant. The expansion is calculated as: 

d 



d = u a a = ±w(r)- 3/2 — (w(r)^/E 2 - w{r)) (18) 



6 Notice that only in Schwarzschild spacetime (which has the property gtt9rr = — 1 )the choose of E = 1 is 
a good coordinate's representation and in general , specially in the form of a cylindrically symmetric LB 
metric, we must keep E ^ 1 to avoiding of occurrence an unphysical radial velocity 
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Not surprisingly, the congruence is diverging (6 > 0) if the geodesies are outgoing, and 
converging (6 < 0) if the geodesies are ingoing.The rate of change of the expansion is 
calculated as: 

d6 _ d6 dr _ ^ r 
dr dr dr 

and the result is: 

fr = { W)~ 1) ^ [ " (rr3/2 ^ ( " (r) V / ^ 2 - W(r))] (19) 
Substituting w(r) and performing differentiation , in terms of a new variables 

w ( r ) / 3/77^ n \ I LI 1 

x = e (V46 2 , 1), \b\ < - 
E 2 2 

E 

y = V x 6 — Ab 2 x 3 
"J 



We have: 

dB_ = A $ b (x,y) 
dr 2 x(l — x) 

Where, 



(20) 



.. . , (27 x 2 - Ahx + 20)y - 366V - 46a; 4 + 27 x b + Ub 2 x + 20x 3 - 326 2 . , 
^ = 3(^ (21) 

This polynomial has at most two real roots, as can be seen from the fact that its second 
derivative is always positive and it vanishes whenever there is a double root. As mentioned 
before, in order to ensure reality, the root point must be at the right of the outer region 
w(r) < l.For |6| < 1/2, < x < 1 this function always poss negative values. For ^ = this 
entails to look for an stationary point of the equation ^(x,y) = 0. This equation can not 
be solved. But for suitable values of b this polynomial in x has a second derivative that is 
everywhere positive regardless of the value of b, implying that it may have at most two real 
roots that can be identified with two event horizons.li one take 6 = 0, when the minus sign 
is chosen, the metric (15) corresponding to the LCA case already analyzed in the previous 
work[l]. These choose will coincide when the roots become w(r) = 0.377, 1.178 and a single 
double root w(r) = (which is not accessible), by w(r) — 1.178 we have 

r = 0.0273a 



For a congruence of null geodesies which must be hypersurface orthogonal (in the manner 
that is discussed in part (a) ) and (NEC) hold, we must have ^| < where 9 = k a a where k a 
is the tangent vector field. 

For dip = dz = the (LB) line element reduces to: 
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ds 2 = -w(r)(dt - - r =){dt + - 1 ==) 
y/w(r) y/w{r) 

The displacements will be null if ds 2 = 0. We define two null coordinates , 

u = t — r 
v = t + r* 

which as usual we introduced a tortoise coordinate 7 

f dr L 1 1 7 2 

r — — = ae a r (-,-,-, — c e ° j 
7 V^H V 6'3'6' ; 

Easily we find that on out-going null geodesies w = constant and similarly t> = constant 
on in-going ones. The following vectors are n«//, 

K ut = -d a u 

K = -9aV 

They both satisfies the geodesic equation with +r as an afline parameter for k^ ut and — r 
for k™. The congruences are clearly hypersurface orthogonal. Expansion(s) are calculated: 

d6 1 rr: 7~ ... . 3w'(r) 2 n 

— = — -a/£ 2 - w r )[w" r) - -—7^- 22 

This function never vanishes and remains always negative. We can construct a similar func- 
tion as(15) for it and determining the sign of it. Another simple method is drawing a graph 
for Appling any of this two methods prove that this function is negative every where. 

SUMMARY 

In this short report we found the unique exact solution for field equations containing a 
massless scalar field and a cosmological constant term . We checked the stability ,also we 
showed that there is some restriction on the energy conditions. This exact solution satisfy 
energy conditions and the validity of focusing theorem proved directly. 
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